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ABSTRACT. 



The chimera overset grid method is reviewed and discussed relative to other procedures for simulat- 
ing flow about complex configurations. It is argued that while more refinement of the technique is 
needed, current schemes are competitive to unstructured grid schemes and should ultimately prove 
more useful. 


INTRODUCTION AND OVERVIEW 


There are currently two mainstream approaches for computing flow fields in which geometry 
imposes complex boundary conditons - composite structured grid schemes and unstructured grid 
schemes. In my assessment of the literature, unstructured grid methods are generally considered 
to be more versatile and easier to adapt to complex geometry while composite structured grid 
methods are generally considered to use more efficient numerical algorithms and require less com- 
puter memory. But with either composite-structured or unstructured grids, the capability to solve 
flow about complex configurations has been suitably demonstrated. Both pure-strain approaches 
have their strengths and weaknesses. Hybrid schemes which incorporate the best features of both 
have already appeared[l-3] and will likely become more important in treating flow about complex 
geometries. 

The chimera[4] and similar methods which use overset grids [5-29] are generally classed into 
the composite structured grid catagory, because these approaches clearly grew out of an attempt 
to generalize body conforming structured grid schemes to treat more complex situations. The 
chimera approach uses a composite of overset structured grids to resolve geometry, flow features, 
or permit more efficient flow solvers. While chimera generally employs composite structured grids, 
the connectivity of the overset structured grids is itself unstructured. 

The chimera approach has been used to compute inviscid and high Reynolds flow about com- 
plex configurations^. f. [7, 8, 14-19, 24-26]), and it has even been demonstrated for unsteady three 
dimensional viscous flow problems in which one body moves with respect to another [26]. The via- 
bility of this approach is perhaps best illustrated by the fact that this progress has been made by a 
relatively small group of researchers. Nevertheless, chimera is sometimes viewed as an intermediate 
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solution approach, one which will ultimately be replaced by the unstructured grid method. With 
further examination, however, the chimera approach may be found to have more versatility than 
current unstructured schemes because, while grids can be abutted together like patches, they can 
also be overset. Oversetting can be somewhat foreign to finite volume and finite element methods, 
but oversetting can be useful. In overset schemes, intermediate boundary curves can be placed 
arbitrarilly. In overset regions, the possibility also exists of impressing solutions from one domain 
onto another via forcing functions rather than only through boundary interfaces. These features 
can and have been exploited in several ways. 

Overset grids allow structured grids to be used without excessive distortion or inefficient use 
of grid density. Consequently, efficient numerical methods can be used which depend on structured 
grids such as alternating direction implicit schemes and parabolized Navier Stokes procedures. 
Moreover, numerical schemes that use structured-grids generally require less computer storage and 
are better suited to vectorizated computers. Arbitrary placement of intermediate boundaries can 
greatly simplify the task of structured grid generation. It allows, for example, the use of hyperbolic 
grid generation procedures[30,31] which do not conform to boundary value constraints but which 
generate nearly orthogonal grids with excellent mesh spacing control. Overset structured grids 
have also been used as a solution adaption procedure[21,22], and overset grids can be positioned 
in the field simply to implement a special solver, most of which require some kind of coordinate 
alignment, (usually to streamlines in flow field simulation). Overset grids allow one body to move 
with respect to another without regriding [10, 26] or placement of new bodies into the domain without 
regenerating the entire mesh. 

As noted earlier, the chimera method is an outgrowth of trying to generalize a powerful solution 
approach, the body conforming structured grid method, to more complex situations. The method 
is proven, but far from mature. There are weaknesses which must be removed if chimera is to 
remain competitve with unstructured grids. 

There are two main criticisms leveled against the current implementations of the chimera 
method. One is the bookkeeping-like complexities associated with connecting overset grids together. 
In fact, the bookkeeping with a chimera scheme is similar to that associated with an unstructured 
grid - easier because the connectivities can be made using structured grid data but more difficult 
because of multiple oversetting. However, only a few reseachers have worked on this problem as 
opposed to the myriad that have worked on unstructured grids. Consequently, there is is less 
available software for overset structured grids. 

The other criticism is that in most simulations of complex flows the solutions on the overset 
grids are merged using simple interpolation. The fact that interpolation is generally used to connect 
grids implies that conservation is not strictly enforced. For most practical applications it is difficult 
to devise a situation where this error is overall-significant since conservation is strictly maintained 
at all points in the domain (assuming the solver is a conservative one) except at a small number 
of interface boundaries. Nevertheless, this source of error has to be eliminated, and conservative 
interface schemes for overset grids have been devised[23,27]. Refinement and simplification of these 
techniques are still warranted for use with three dimensional flow solvers, however. 

In the remainder of this paper, a brief review of the chimera scheme is given, and a few results 
from previously presented space shuttle flow field simulations are used to indicate current status. 
Some possible future directions for the chimera scheme are then indicted, followed by concluding- 
remarks. 
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BACKGROUND AND PROGRESS 


The chimera composite grid discretization method is a domain decomposition approach which 
uses overset body-conforming grids. In this approach, each component of a configuration is grided 
separately and is overset onto a major grid to form the complete model. The major grid is stretched 
over the entire field, and is often generated about a dominant boundary or body surface. Minor 
grids are used to resolve features of the geometry or flow that are not adequately resolved by the 
major grid, and are overset on the major grid without requiring mesh boundaries to join regularly. 

For example, Fig. 1 shows surface grids generated for the integrated space shuttle configu- 
ration in its ascent mode. The configuration shown has simplified attach hardware, and various 
protuberances such as the external fuel lines and even the orbiter vertical tail have been neglected. 
A grid is then independently generated about each component. A composite grid is then formed by 
superimposing all grids together. The body-conforming grids used for each component are shown 
in Fig. 2 at their respective planes of symmetry. Here the external tank (ET) grid is treated as 
the major grid and is extended to the far field. Figure 3 shows nearby £ = constant planes for the 
orbiter and ET projected onto an x=constant plane. Whenever points of a grid, say grid 1, fall 
within the body-boundary of another grid, say grid 2, the points of grid 1 are cut out forming a 
hole in grid 1. The hole-boundary data of grid 1 are then supplied from grid 2. Hole grid points 
have been removed from view in Fig. 3. 

Software to interconnect the grids is needed to ascertain when points of one grid fall within a 
body boundary of another (grid hole points) and to supply pointers so that one grid can provide 
boundary data to another. Various algorithms have been devised for performing these tasks auto- 
matically [4,8,14-16,20,29]. For the illustrated shuttle grids, the code Pegasus [8,14,16,29] (provided 
and maintained by CALSPAN of AEDC) has been used. General software for this problem has also 
been developed in [20] and includes an interactive workstation demonstrator for two dimensional 
grids. 

A flow simulation code developed for a single general curvilinear grid can be readily adapted 
for composite overset grids. One simply sets flags to blank out hole points and supplies a control 
program that calls in grids and interface routines. For example, the structured- grid, implicit, 
approximately-factored F3D scheme [32,33] for the thin-layer Navier-Stokes equations 

d T Q + d^F + d v G + d c H = Re- 1 d c S (1) 

was easily modified for chimera overset grids as 


[/ + i b hSl(A + ) n + i 6 M c C n - i b hRe~ 1 S < J~ x M n J - i b hDi \ c 

x [/ + i b hS( ( A~) n + i b h5 v B n - i b hDi\ n ] A Q n = (2) 

- i b At{Sl(F+) n + 6%(F-) n + S v G n + - Re~ l 6 Q S n + (D e \ v + D e \ c )Q n } 

Here introduction of the flag i b accommodates the possibility of having arbitrary holes in the grid. 
(The hole includes hole-boundary points which are later updated by interpolating the solution 
from the overset grid which created the hole.) The array of values i b is defined such that i b = 1 
at normal grid points and i b = 0 at hole points. Thus, in Eq. (2) when i b = 1 the normal scheme 
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Fig. 1 Simplified configuration and surface grid point distributions for the integrated space 
shuttle. 




is maintained, but when i b = 0 the scheme reduces to A Q n = 0 or Q n+1 = Q n and thus Q is 
not changed at a hole point. For the most part i b is coupled to the time step (h = At or (Ai)/2) 
and is trivial to implement into the coding. (Difference operators that use more than one point to 
either side require some additional coding modification, see Ref. [24]) By using the i b array it is 
not necessary to provide branching logic to avoid hole points, and computer vectorization is not 
inhibited. 

The F3D algorithm was implemented on composite grids by writing a control program which 
at each “time step” fetched a grid and its data from an isolated large memory into a working 
memory. Boundary interface arrays that store grid interconnect data, Qbc > are also fetched. The 
Qbc array holds overset-grid boundary values for the current grid which are supplied from the 
other grids and is a relatively small array. Because the hole boundaries are arbitrarily located, the 
Qbc array has pointers much like those used with unstructured grids. The solution on the current 
grid is then updated or advanced in time. Overset boundary data that the current grid sends to 
other grids are then found by interpolation and loaded into QbCi and all arrays are sent back to 
the isolated large memory. The next grid is then fetched, and so on. 

To illustrate this capability, calculated results for the integrated space shuttle vehicle are taken 
from Refs. [24,25] and are reproduced in Figs. 4 and 5. These figures show comparisons between 
computational and experimental data for = 1.05 at an angle of attack, a = —3°, and using 
the wind tunnel Reynolds number Re = 4.0 x 1 0 6 /ft for the computations. Shaded surface pressure 
coefficient comparisons between the computation and wind tunnel data [34] are shown by Fig. 4. 
This kind of comparison is possible because the 3% scale wind tunnel model was instrumented with 
1538 pressure taps. Mach contours in planes of symmetry of the ET and solid rocket booster (SRB) 
are also shown in Fig. 4, and are used to highlight the SRB plume which was modeled as a hot-air 
jet. A limited amount of flight test data [35] are also available for comparison, and Fig. 5 shows 
pressure comparisons between computation, flight, and wind tunnel data taken along the side of 
the fuselage. This computation required about 15 hours of computer time using a single processor 
of the CRAY 2, and employed a composite grid containing one million points distributed over seven 
distinct grids. Additional details of the experimental comparisons (and some disclaimers) are given 
in Refs. [24,25]. 

With the composite overset grid approach it is possible to move one body with respect to 
another without regriding at each time step advance of the flow field. Demonstration of this 
capability for SRB staging from the shuttle is presented in Ref. [26] . 


FUTURE DIRECTIONS 


As noted previously, two main criticisms can be leveled against the chimera approach: 1) the 
complexity of the interconnectivity is perhaps as difficult as dealing with an unstructured grid, 
and 2) nonconservative interpolations to update interface boundaries are often used in practical 
three dimensional computations. To indicate the complexity of the interconnectivity, it should be 
remarked that some of the attach-hardware used in the space shuttle simulations is not actually 
attached. The attach-hardware was floated between the body elements. This is because the algo- 
rithms devised to impose grid connectivity are not accurate enough when the refined grids used for 
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b) orbiter gnd 


Fig. 3 Cross-section of grids showing holes, a) ET grid, blOrbiter grid 
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Fig. 4 Comparison of pressure coefficient between computation (top) and wind tunnel (bottom) 
Moo = 1.05, a — —3°, and Re = 4.0 x 10 6 //f (3% model). 
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Fig. 5 Comparison of C p from computation (-), wind tunnel (o), and flight test (V right side, A left 
side) along the (f> = 70° line of the orbiter fuselage, = 1.05 and a — -3°, Re = 4.0 x 10 6 //t, 
and 10° /9° elevon deflection. 



high Reynolds number viscous flow simulation are encountered. Because the normal grid spacing 
is so fine, small errors in interpolation for the body surface can cause the test which identifies 
when points lie inside a another body surface to fail. Specifically, two viscous boundary layer grids 
emanating off the same surface cannot fall within that surface, yet due to the fineness of the grids, 
a viscous flow field point above the surface of one grid may be judged to be inside the same body 
surface of another grid because of interpolation error. Special logic can be used to exclude this case, 
and the problem is not encountered for inviscid simulations. There are a variety of ways to deal 
with this problem, they include more consistent interpolation schemes, introducing special viscous 
surface grids, alternate tests, and so on, but current software is inadequate and needs refinement. 

Several approaches have been proposed to treat the problem of nonconservative interfaces, 
but they have not been implemented into the F3D code. It should be remarked that while the 
space shuttle simulations have used simple interpolation procedures because of their robustness, 
CALSPAN simulations for the last several years[16-18,28,29] have implemented an unpublished 
idea of Benek and use interpolants of delta quantities, specifically, Q n+1 — Q n . Interpolating this 
quantity on interface boundaries ensures space-time conservation over the global field, but the 
ultility (or penalty) of this approach has not yet been rigourously examined. 

A fairly obvious way to ensure interface conservation would be to introduce an unstructured 
flow solver in the vicinity of the interface boundaries. Already in chimera, primitive elements of 
an unstructured grid solver exist in the form of pointers and grid interfacing arrays which transfer 
interpolated values of the solution from one grid to the next. Some care would have to be taken, but 
an explicit differencing of the governing equations using unstructured data could substitute for the 
interpolation process. The chimera would then mimick a hybrid structure-unstructured approach 
much as in Refs. [1-3]. 

Finally, the chimera framework lends itself to the construction of a general-purpose flow code 
that can optionally take advantage of approximate solution methods, and some preliminary work 
has been carried out in this area. For example, the F3D flow simulation code used for the previously 
described shuttle work already has options (at various levels of maturity) to use either explicit or 
implicit solution algorithms as well as a semimarching scheme for predominately steady supersonic 
flow. In addition, a fortified option [36] of the basic algorithm is available. To support this option 
and to provide diagnostics of Navier-Stokes solution accuracy, a three-dimensional boundary-layer 
code in arbitrary general coordinates has also been included. 


CONCLUDING REMARKS 


Overset grid schemes such as chimera grew out of efforts to adapt body conforming structured 
grid methods to more complex boundaries. Although these schemes have not received the attention 
that unstructured (or patched) grid methods have received, they have proven to be competitive, 
and are likely to prove to be more powerful. Incorporation of unstructured grids into the overall 
chimera framework is quite feasible, and offers a relatively low risk route to a hybrid structured- 
unstructured simulation code as well as a fall back position in the (unlikely) event that chimera 
does not live up to its expectations. 
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